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Abstract 


Association  of  radar  and  ESM  tracks  is  an  important  component  for  multisensor- 
multitarget  tracking  systems.  It  is  essential  to  the  overall  performance  of  a  track  fusion 
process.  In  this  report,  we  present  an  ESM/track  association  algorithm  in  a  new  modified 
polar  coordinate  (MPC)  system.  The  algorithm  is  based  on  a  theoretical  performance  analysis 
of  ESM/radar  track  association  techniques  in  the  MPC.  Formulas  for  computing  the  correct 
and  false  association  probabilities  with  and  without  state  estimation  biases  are  derived.  We 
also  discuss  the  effect  of  the  number  of  observable  states  on  association  performance.  The 
proposed  algorithm  uses  a  maximum  likelihood  (ML)  algorithm  for  estimating  the  target 
states,  of  which  the  observable  state  estimates  are  subsequently  used  for  the  ESM/radar  track 
association.  Computer  simulations  are  used  to  demonstrate  the  performance  and  improvement 
of  the  proposed  ESM/radar  track  association  method. 


Resume 

L’association  de  pistes  radar  et  SME  est  une  composante  importante  des  systemes 
multicapteurs  de  poursuite  de  cibles  multiples.  Cette  association  est  essentielle  au  rendement 
global  du  processus  de  fusion  de  pistes.  Dans  le  present  document,  nous  presentons  un 
algorithme  dissociation  de  pistes  SME/radar  dans  un  nouveau  systeme  de  coordonnees 
polaires  modifie  (MPC).  L’algorithme  est  base  sur  une  analyse  de  rendement  theorique  des 
techniques  dissociation  de  pistes  SME/radar  dans  le  systeme  MPC.  Les  formules  servant  a 
calculer  les  probability  dissociation  correctes  et  erronees  avec  et  sans  les  biais  d’estimation 
d’etat  sont  derivees.  Nous  discutons  egalement  des  effets  d’un  bon  nombre  d’etats 
observables  sur  le  rendement  de  T association.  L ’algorithme  propose  utilise  un  algorithme  de 
maximum  de  vraisemblance  (MV)  pour  estimer  les  etats  des  cibles.  Ces  etats  observables 
estimes  sont  ensuite  utilises  pour  Tassociation  de  pistes  SME/radar.  Les  simulations 
informatiques  servent  a  demontrer  le  rendement  et  1’ amelioration  de  la  methode  dissociation 
de  pistes  SME/radar. 
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Executive  summary 


This  report  focuses  on  the  integration  of  radar  and  electronic  support  measures  (ESM) 
sensors,  in  particular,  the  association  of  ESM/radar  tracks.  The  work  is  motivated  by  the 
requirement  of  the  DRDC  TDP  SISWS  (Shipbome  Integration  of  Sensor  &  Weapon  System) 
project,  and  may  contribute  to  the  sensor  manager  part  of  the  project  in  the  future. 

Fusion  of  target  tracks  from  multiple  similar  or  dissimilar  sensors  is  an  effective  method  to 
improve  tracking  accuracy  and  target  acquisition  rates.  ESM  sensor  and  radar  are  the  two 
most  popular  types  of  sensors  that  have  been  deployed  in  air,  sea  and  land  surveillance 
systems.  The  information  obtained  by  radar  and  ESM  sensors  is  complementary,  and 
integrating  radar  and  ESM  sensor  outputs  will  help  to  increase  the  likelihood  of  target 
acquisition  and  to  improve  the  tracking  quality. 

Track  association  is  a  pre-requisite  for  ESM/radar  fusion  and  is  essential  to  the  overall 
performance  of  tracking  systems.  One  of  the  uniqueness  of  ESM/radar  track  association  is  that 
the  only  common  parameter  between  the  ESM  sensors  and  radar  is  the  bearing  measurements. 
In  addition,  the  bearing  measurements  are  not  synchronized.  Usually  bearing-only  tracking  is 
required  for  ESM  sensors  for  track  association.  However,  bearing-only  tracking  usually 
suffers  from  the  observability  problem,  i.e.,  the  target  states  may  not  be  observable  when  only 
bearing  measurement  is  available.  In  order  to  make  the  all  of  the  target  states  observable,  it 
requires  that  the  ESM  sensor  have  a  motion  with  non-zero  derivative  of  higher  order  than  the 
target,  i.e.,  the  platform  has  to  maneuver  more  rapidly  that  the  targets.  This  is  a  rare  scenario 
in  practice,  especially  for  shipborne  ESM  sensors.  Non-observable  states  will  result  in 
inaccurate  and  biased  estimates,  which  in  turn,  will  severely  deteriorate  the  performance  of 
ESM/radar  track  association. 

In  this  report,  we  present  ESM/track  association  algorithm  in  a  new  modified  polar 
coordinate  (MPC)  system  based  on  a  theoretical  performance  analysis  of  ESM/radar  track 
association  techniques  in  the  MPC.  A  conclusion  from  the  performance  analysis  is  that ,  under 
the  Gaussian  assumption,  using  more  observable  states  would  be  able  to  improve  the 
performance  of  ESM/radar  association.  In  the  proposed  algorithm,  a  state  vector  is  used  which 
includes  extra  observable  states  in  addition  to  the  first  three  states  in  Aidala’s  MPC.  The 
bearing-only  tracking  of  the  ESM  sensors  is  implemented  using  a  maximum  likelihood  (ML) 
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algorithm.  Computer  simulations  show  that  the  association  performance  can  be  improved 
effectively  by  using  more  observable  states  in  the  new  MPC  in  different  scenarios. 


Y.  Zhou,  W.  Li  and  H.  Leung.  2004.  Maximum  Likelihood  Based  ESM/Radar  Track 
Association  Algorithm  in  a  New  Modified  Polar  Coordinate.  DRDC  Ottawa  TM  2004- 
255.  Defence  R&D  Canada  -  Ottawa. 
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Sommaire 


Le  present  document  porte  principalement  sur  1’ integration  de  capteurs  radar  et  de  capteurs 
mesure  de  surveillance  electronique  (MSE),  et  plus  particulierement  sur  l’association  de  pistes 
MSE/radar.  Les  travaux  sont  justifies  par  les  besoins  du  projet  sur  l’lntegration  de  capteurs  et 
de  systemes  d'armes  embarques  (SISWS)  du  Programme  de  demonstration  de  technologies 
(TDP)  de  RDDC,  et  pourront  servir  plus  tard  a  la  partie  du  projet  portant  sur  la  gestion  des 
capteurs. 

La  fusion  de  pistes  de  cibles  provenant  de  plusieurs  capteurs  similaires  ou  differents  est 
une  methode  efficace  pour  augmenter  la  precision  des  poursuites  et  les  taux  d’ acquisition  de 
cibles.  Les  capteurs  SME  et  radar  sont  les  deux  types  de  capteurs  les  plus  utilises,  et  ils  sont 
deployes  dans  les  systemes  de  surveillance  aerienne,  maritime  et  terrestre.  Les  informations 
obtenues  par  les  capteurs  radar  et  SME  sont  complementaires,  et  l’integration  des  sorties  des 
capteurs  radar  et  SME  aidera  a  augmenter  la  vraisemblance  de  V  acquisition  de  cibles  et  la 
qualite  des  poursuites. 

L’association  de  pistes  est  un  prealable  a  la  fusion  de  pistes  SME/radar,  et  elle  est 
essentielle  au  rendement  global  des  systemes  de  poursuite.  Une  des  specificites  de 
l’association  de  pistes  SME/radar  est  que  le  seul  parametre  commun  entre  les  capteurs  SME  et 
radar  est  la  mesure  de  gisement.  De  plus,  les  mesures  de  gisement  ne  sont  pas  synchronises. 
Generalement,  la  poursuite  en  gisement  seul  est  necessaire  pour  les  capteurs  SME  aux  fins 
dissociation  de  pistes.  Toutefois,  la  poursuite  en  gisement  seul  est  generalement  touchee  par 
le  probleme  d’ observability,  c’est-a-dire  que  les  etats  de  cible  peuvent  ne  pas  etre  observables 
lorsqu’on  dispose  seulement  des  mesures  de  gisement.  Afin  que  tous  les  etats  de  cible  soient 
observables,  il  faut  que  le  capteur  SME  presente  un  mouvement  avec  une  derivee  non  nulle 
d’un  ordre  superieur  a  celui  de  la  cible,  c’est-a-dire  que  la  plate-forme  doit  se  deplacer  plus 
rapidement  que  la  cible.  Ce  scenario  est  plutot  rare  en  pratique,  particulierement  pour  les 
capteurs  SME  de  navire.  Des  etats  non  observables  resultent  en  estimations  imprecises  et 
faussees,  qui  a  leur  tour  peuvent  grandement  deteriorer  le  rendement  de  l’association  de  pistes 
SME/radar. 

Dans  le  present  document,  nous  presentons  un  algorithme  dissociation  de  pistes 
SME/radar  dans  un  nouveau  systeme  de  coordonnees  polaires  modifie  (MPC)  base  sur  une 
analyse  de  rendement  theorique  des  techniques  dissociation  de  pistes  SME/radar  dans  le 
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systeme  MPC.  La  conclusion  que  nous  pouvons  tirer  de  l’analyse  de  rendement  est  que,  selon 
l’assomption  de  Gauss,  l’utilisation  d’etats  plus  observables  permettrait  d’ameliorer  le 
rendement  de  l’association  de  pistes  SME/radar.  Dans  Ealgorithme  propose,  on  utilise  un 
vecteur  d’etat  qui  inclut  des  etats  observables  supplementaires  en  plus  des  trois  premiers  etats 
du  systeme  MPC  d’Aidala.  La  poursuite  en  gisement  seul  des  capteurs  SME  est  mise  en 
application  a  l’aide  d’un  algorithme  de  maximum  de  vraisemblance  (MV).  Les  simulations 
informatiques  montrent  que  le  rendement  de  1’  association  peut  etre  ameliore  de  fagon  efficace 
si  on  utilise  plus  d’etats  observables  avec  le  nouveau  systeme  MPC  dans  differents  scenarios. 


Y.  Zhou,  W.  Li  and  H.  Leung.  2004.  Maximum  Likelihood  Based  ESM/Radar  Track 
Association  Algorithm  in  a  New  Modified  Polar  Coordinate.  DRDC  Ottawa  TM  2004- 
255.  R  &  D  pour  la  defense  Canada  -  Ottawa. 
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1.  Introduction 


Fusion  of  target  tracks  from  multiple  similar  or  dissimilar  sensors  is  an  effective  method  to 
improve  tracking  accuracy  and  target  acquisition  rates  [1].  In  practice,  we  are  usually 
interested  in  fusing  radar  and  electronic  support  measures  (ESM)  sensor  outputs  because  of 
two  reasons.  First,  radar  and  ESM  sensors  are  the  two  most  popular  types  of  sensors  that  have 
been  deployed  in  air  surveillance  systems  [2]-[5].  Secondly,  the  outputs  of  radar  and  ESM 
sensors  are  complementary.  Radar  usually  reports  kinematics  parameters  of  the  targets 
(position,  speed,  acceleration  with  covariance  information)  while  ESM  sensors  measure  the 
attributes  of  the  target  emitters,  which  include  information  on  emitter  RF,  PRI,  PW  as  well  as 
their  bearings.  The  information  obtained  by  radar  and  ESM  sensors  is  complementary,  and 
integrating  radar  and  ESM  sensor  outputs  will  help  to  increase  the  likelihood  of  target 
acquisition  and  to  improve  the  tracking  quality.  The  overall  outcome  of  fusion  is  improved 
situational  awareness  and  more  accurate  and  complete  composite  target  track  files. 

In  order  to  fuse  radar  and  ESM  sensor  information,  radar  and  ESM  tracks  need  to  be 
associated,  i.e.,  tracks  from  the  two  sensors  that  represent  the  same  emitter  be  recognized.  It  is 
a  pre-requisite  for  ESM/radar  fusion  and  essential  to  the  overall  performance  of  tracking 
systems.  ESM/radar  track  association  is  unique  in  that  the  only  common  parameter  between 
them  is  bearing  measurement.  In  addition,  the  tracks  from  radar  and  ESM  sensors  are  not 
synchronized  in  time  because  ESM  tracks  mainly  depend  on  the  illumination  of  the  target 
emitters  and  the  scheduling  of  the  receivers.  Since  the  only  positional  information  that  the 
ESM  sensors  observe  is  the  bearing  measurement,  bearing-only  tracking  is  required  for  ESM 
sensors  for  track  association.  The  association  performance  follows  closely  that  the  bearing- 
only  tracking  at  the  ESM  sensors.  One  issue  with  bearing-only  tracking  is  the  observability, 
i.e.,  the  target  states  may  not  be  observable  when  only  bearing  measurement  is  available.  It 
was  found  that  all  the  target  dynamical  states  couldn’t  be  estimated  accurately  when  the 
targets  are  non-maneuvering.  In  recent  years,  many  bearing-only  tracking  algorithms  have 
been  developed  in  the  literature  [6-12]  to  tackle  the  ovservability  problem.  The  proposed 
approaches  usually  implement  the  bearing-only  tracking  in  different  coordinate  systems 
including  the  polar  coordinate  (PC)  [6-8],  the  Cartesian  coordinate  (CC)  and  the  modified 
polar  coordinate  (MPC)  [11].  The  MPC  approach  by  Aidala  [11]  has  been  shown  to  be  able  to 
overcome  some  shortcomings  of  passive  tracking  associated  with  the  CC.  In  Aidala’ s 
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approach,  the  four  states  in  the  MPC  are  bearing  rate,  ratio  of  range  rate  and  range,  bearing, 
and  inverse  range.  He  showed  that  the  bearing  state  is  always  observable,  and  the  bearing  rate 
and  the  ratio  of  range  and  range  rate  are  observable  almost  everywhere  except  in  some  special 
cases.  However,  to  make  the  inverse  range  observable,  the  sensor  platforms  or  the  targets 
must  be  maneuvering  [13 [-[16].  In  [12],  the  performances  of  some  track  association  methods 
based  on  bearing-only  tracking  in  the  CC  and  the  MPC  were  compared  using  the  computer 
simulations.  In  [9] [17],  theoretical  performance  was  carried  out  for  a  multisensor  track-to- 
track  correlation  technique  in  the  CC,  which  includes  the  correct  and  false  association 
probabilities. 

In  this  report,  we  present  an  ESM/track  association  algorithm  in  a  new  modified  polar 
coordinate  (MPC)  system.  The  algorithm  is  based  on  a  theoretical  performance  analysis  of 
ESM/radar  track  association  techniques  in  the  MPC.  Note  that  the  above-mentioned 
performance  analysis  is  limited  to  particular  state  models  used  in  the  MPC  and  applicable  to 
state  vectors  of  fixed  numbers  of  states.  We  develop  a  more  general  performance  analysis, 
which  is  suitable  for  any  number  of  states.  In  the  analysis,  correct  and  false  association 
probabilities  for  maneuvering  and  non-maneuvering  targets  are  derived.  From  the  analysis,  we 
show  that  the  approach  by  Aidala  always  results  in  biased  estimation  of  the  fourth  state,  i.e., 
the  inverse  range,  due  to  the  non-observability  for  non-maneuvering  targets  and  observers. 
The  effect  of  the  estimation  bias  of  the  inverse  range  state  on  the  association  performance  is 
discussed.  As  a  conclusion,  we  show  that,  under  the  Gaussian  assumption,  using  more 
observable  states  can  be  used  to  improve  ESM/radar  association  performance.  To  counter  the 
observability  of  problem  associated  with  non-maneuvering  targets,  we  propose  an  approach  in 
which  the  state  vector  includes  extra  observable  states  in  addition  to  the  first  three  states  in 
Aidala’ s  MPC.  The  state  vector  has  five  states,  namely,  bearing,  ratio  of  range  rate  and  range, 
bearing  rate,  heading  direction,  and  ratio  of  target  speed  and  range.  The  heading  direction  and 
the  ratio  of  the  target  speed  and  range  were  first  defined  in  [18]  [19]  for  linearly  uniformly 
moving  targets,  and  were  shown  to  be  observable  almost  everywhere  except  for  an  extreme 
condition.  The  bearing-only  tracking  of  the  ESM  sensors  is  implemented  using  a  maximum 
likelihood  (ML)  algorithm.  The  ESM/radar  track  association  is  then  to  convert  the  radar  tracks 
to  the  new  MPC  and  associate  them  with  the  ESM  tracks.  The  CRB  of  the  ML  estimator  is 
also  derived.  It  should  be  mentioned  that,  when  the  targets  or  the  observers  are  maneuvering, 
the  additional  observable  states  can  still  be  used  to  improve  the  track  association  performance. 
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This  report  is  organized  as  follows.  In  Section  2,  we  analyze  the  theoretical  performance 
of  classical  ESM/radar  track  association  techniques  for  maneuvering  and  non-maneuvering 
targets.  The  effect  of  the  number  of  observable  states  on  the  association  performance  is  also 
discussed.  In  Section  3,  the  system  model  for  track  association  in  the  new  MPC  is  introduced. 
Section  4  is  devoted  to  the  ML  estimation  method  for  ESM/radar  sensor  tracking  and 
association  algorithm  proposed  in  the  new  MPC.  Observability  analysis  for  ESM  bearing-only 
tracking  in  the  new  coordinate  system  is  also  given  in  this  section.  Finally,  in  Section  5, 
computer  simulations  are  used  to  demonstrate  the  performance  of  the  proposed  track 
association  algorithm. 
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2.  Performance  of  Calssical  ESM/Radar  Track 
Association Techniques 

In  this  section,  we  discuss  the  association  performance  for  maneuvering  and  non-maneuvering 
targets.  In  particular,  the  effect  of  the  number  of  observable  states  on  association  performance 
is  discussed. 

2.1  Association  performance  for  manoeuvring  and  non¬ 
manoeuvring  targets 

Consider  the  near-constant  velocity  model.  Using  the  Cartesian  coordinate,  the  target  state 
vector  is  given  by  xcc  =  [x,  x,  y,  y]T ,  where  x  and  y  are  the  positions,  and  x  and  y  are  the 
velocity  along  x-  and  y-axis,  respectively.  The  MPC  proposed  by  Aidala  [11]  is  given  by 
xmpc  =  [9,r  /  r,9,\/ r]T  ,  where  9 ,  f  /  r  ,  9 ,  and  1/r  represent  the  bearing  rate,  ratio  of  range 
rate  and  range,  bearing,  and  inverse  range,  respectively.  Assume  that  the  target  state  estimates 
by  radar  and  the  ESM  sensor  are  independent  and  Gaussian  distributed  with  zero  mean  and 
variances  given  by 


^  radar  diag 

^ esm  ~  diag 


o2  c*2  c*2  o2  I 

0 6, radar  ’  °r!r, radar  ’  ° 6, radar  ’  °i  I  r, radar  I 

^6, ESM  5  ^r/r,ESM  >  d0  ESM  ,  S\/r^ESM  J  ? 


respectively.  Let  us  denote  xradar  and  xESM  as  the  state  estimates  by  radar  and  ESM  sensor, 
respectively,  the  association  is  to  determine  whether  they  are  from  the  same  target  or  not.  A 
statistical  test  for  a  four  state  vector  can  be  written  as 


Y~(X radar  XESM  )  radar  +^ESM  )  ( X radar  XESM  )  ‘ 


(1) 


The  analysis  can  be  categorized  into  four  cases:  (1)  two  tracks  are  from  the  same  target  and 
all  target  states  are  estimated  unbiased;  (2)  two  tracks  are  from  the  same  target  and  one  of  the 
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states  estimates  is  biased;  (3)  two  tracks  are  from  different  targets  and  all  target  state  estimates 
are  unbiased;  (4)  two  tracks  are  from  different  targets  and  one  of  the  states  is  biased. 

In  the  first  case,  the  statistic  Y  can  be  shown  to  be  chi-square-distributed  with  n  degrees  of 
freedom  [20] 


Pr(y)  =  - 


1 


2n,1T{n!2) 


y 


«/2-l  -y!2 


y>0, 


(2) 


where  n  is  the  number  of  target  states  used  for  track  association,  and  T(g)  is  the  gamma 
function 


r (q)  =  ^t9~le~,dt,  q>  0, 

r(<?)  =  (<7_l)!>  <7  an  integer,  q  >  0,  (3) 

r(0.5)  =  yjn,  T(1 .5)  =  0.5>/jr. 


Let  the  test  threshold  be  T  .  The  correct  and  false  association  probability  can  be  written  as 


Pc 


\lpY(y)dy  =  J( 


1 


2nllT(n/2) 


■y 


n/2~le~y/2d 


ly, 


(4) 


and 


pf  =J"  pY(y)dy=l~Pc’ 


(5) 


respectively.  For  even  n,  the  correct  association  probability  (4)  can  be  simplified  as 


P=\-  e 


-r/2”^-1  1 


k=o  k\ 


{  t\ 


(6) 


In  the  second  case,  assume  that  the  /?th  state  is  unobservable  with  an  estimation  bias 
denoted  by  A mn .  It  is  also  assume  that  the  rest  (n- 1)  estimates  are  independently  Gaussian 
distributed  with  zero  mean.  It  can  be  verified  that  the  Y  follows  a  non-central  chi-square 
distribution  with  n  degrees  of  freedom 
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pY(y)  = 


f  \(«-2)/4 

7 


\AmnJ 


(J,+AW',)/2/„/2-l(^Aw„),7>0 


(V) 


where  Ia  (x)  denotes  the  a  -th  order  modified  Bessel  function  defined  by 


\a+2k 


ia(x)=t  (x/2)° - ,  x>0. 

“  i^ok\Y(a  +  k  +  l) 


(8) 


The  correct  association  probability  can  be  obtained  as 


Pc  =  ft Pr(y)<ty  =  lo  (V2) 


f  \(n- 2)/4 

y 

\Amn  J 


-(W)/2/„/2_l(^Am„)^, 


(9) 


and  the  false  association  probability  can  be  computed  by  (5). 

In  the  third  case,  consider  the  normalized  element  (xradar  i  -  xESM  i)2  / [52adard  +  82ESM  i  j  .  It 

can  be  verified  that,  under  the  Gaussian  assumption,  the  normalized  element  has  a  variance  of 
unity.  Assuming  that  the  normalized  element  has  a  mean  denoted  by  mi ,  the  statistics  Y  can 
be  verified  to  follow  a  non-central  chi-square  distribution  with  n  degrees  of  freedom 


Pr(y): 


f  \(n-2)/ 4 

_y_ 

4n  J 


^y+qlV2Kn-A^yqn\ 


(10) 


where 


q2n  -Tmf  . 

i= 1 


The  correct  and  the  false  association  probability  are  given  by 


(ii) 
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0-2)/4 


(12) 


Pc  =  \rPY^y)dy  =  J7(1/2) 


0-(y+qn)12 


n!  2-1 


(\fyqn)dy  > 


and 


/V  =  f  PY(y)dy=l~Pc » 


(13) 


respectively. 

In  the  fourth  case,  the  tracks  are  from  different  targets  and  the  estimate  of  the  /?th  state  has  a 
bias  Amn  .  The  correct  and  the  false  association  probability  can  be  computed  by  (12)  and  (13), 

respectively,  with  qn  given  by 


n- 1 

< In  =2>;  +K  + 

1=1 


(14) 


2.2  Effect  of  number  of  observable  states  on  association 
performance 


When  a  target  state  is  unobservable,  its  estimate  will  be  biased.  A  biased  state  estimate  will, 
in  turn,  deteriorate  the  association  performance  in  terms  of  correct  association  probabilities. 
Most  association  algorithms  have  been  based  on  the  use  of  the  observable  state.  Obviously, 
the  use  of  different  coordinate  systems  has  an  impact  on  the  performance  of  association.  In  the 
following,  we  discuss  how  the  number  of  observable  states  would  affect  the  performance  of 
association. 

Consider  different  sets  of  observable  states  for  ESM/radar  track  association.  When  only  the 
bearing  is  used  for  association,  the  statistical  test  variable  Y  is  chi-square  distributed  with  a 
degree  of  freedom  of  one.  Let  Tx  denote  the  testing  threshold.  When  more  observable  states 
are  used,  the  new  test  variable  can  be  written  as 


Y  —  (  *  radar  *ESM  )  radar  +  ^ESM  )  (  A 


radar 


VESM 


). 


(15) 
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which  can  be  verified  to  be  a  chi-square  random  variable  with  zero  mean  and  variance  l/n.  If 
the  same  test  threshold  7J  is  used,  the  correct  association  probability,  assuming  that  the  tracks 
are  from  the  same  target,  is  given  by 


Pc =S  PY^y)dy=\l' 


1 


(l/«)"2"/2r(n/2) 


rJ 


nl2-\e-ynl2d. 


y- 


(16) 


For  an  integer  l  =  nil,  (16)  can  be  written  in  a  simpler  form  as 


t!  1 


pc=l-e-t'n"£- 
k=o  k\ 


f  7 


v  2  J 


(17) 


From  (17),  it  can  be  seen  that,  when  the  number  n  of  the  observable  states  increases,  the 
correct  association  probability  improves. 

When  the  tracks  are  from  different  targets,  it  can  be  verified  that  7  is  a  non-central  chi- 
square  variable  with  a  nonzero  mean  and  variance  1  In.  The  false  association  probability  is 
given  by 


Pf  =  \l'pY(y)dy  =  \l 


\{n~2)!A  /  2  \ 

^  e  qn+yn/2ln/2-l(^n<ln)dy’ 


Wn) 


(18) 


which  depends  on  both  the  number  of  observable  states  n  and  qn .  The  parameter  qn  is 
determined  by  the  ratio  of  the  state  difference  and  the  state  estimation  variance.  The 
probability  difference  of  falsely  associating  tracks  from  different  targets  with  (ft+1)  and  n 
observable  states  can  be  expressed  as 


Apf=C 


n  +  \ 

(  V 

y 

2 

Kq2n+  m2n+x , 

\(«-l)/4 


-(ql+ml+x+y\n+\)l2  r—  /“T 
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dy, 
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where  mn+l  represents  the  ratio  of  the  (n+  l)th  observable  state  difference  for  two  tracks  from 
different  targets  and  the  (?2+l)th  state  estimation  variance.  There  is  no  explicit  formula  for 
determining  the  value  of  mn+l  that  satisfies  A pf  <  0 .  If  mn+l  is  large  enough  (as  will  be 

shown  in  the  next  section),  the  false  association  probability  will  decrease.  Another  criterion  is 
defined  here  to  describe  the  averaged  correct  association  probability  of  associating  tracks 
from  the  same  target  and  from  different  targets, 

PraSe=(Pc+0-Pf))/2.  (20) 

To  ensure  that  pc  increases  and  pf  decreases  at  the  same  time  when  using  one  more 
observable  state  (i.e.,  the  (n+  l)th  state),  the  ratio  mn+l  of  the  («+ 1)  state  difference  and 
estimation  variance  must  satisfy  A pf  <  0 . 

2.3  Simulation  results 

In  this  section,  we  use  computer  simulations  to  demonstrate  the  performance  of  classical 
ESM/radar  track  association  techniques.  Figure  1  shows  the  correct  association  probabilities 
when  bearing/inverse  range  or  the  four  states  in  Aidala’s  MPC  under  different  cases.  In  Figure 
1(a),  the  radar  and  ESM  tracks  are  originated  from  the  same  target  while  in  Figure  1(b)  they 
are  from  different  targets.  When  the  observer  is  made  to  be  optimally  maneuvering  [16],  the 
four  states  are  observable  for  bearing-only  tracking  by  the  ESM  sensor.  If  the  target  and  ESM 
sensor  are  not  maneuvering,  the  inverse  range  is  not  observable  and  will  result  in  an 
estimation  bias.  Figures  1(a)  and  (b)  show  the  correct  association  probability  pc  in  the 
absence  of  bias  for  the  inverse  range,  w  =  0,  when  the  tracks  are  from  same/different  targets, 
respectively.  It  is  shown  that  the  correct  association  probabilities  deteriorate  with  the  increase 
of  the  estimation  bias  w  of  the  inverse  range.  Therefore,  the  unobservable  states  are  not 
preferred  for  track  association  in  the  classical  track  association  techniques. 

Figure  2  shows  the  effect  of  the  number  of  observable  states  on  the  performance  of  the 
ESM/radar  track  association  techniques.  The  testing  threshold  is  fixed  at  T  =  3.84.  From 
Figure  2(a),  we  can  observe  that,  when  the  tracks  are  from  the  same  target,  the  correct 
association  probability  increases  when  extra  observable  states  are  used.  Figure  2(b)  shows  the 
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false  association  probability  versus  the  ratio  of  the  (/?+ 1  )th  state  difference  and  estimation 
variance  (n  =  1  here).  The  ratio  of  the  first  state  difference  and  estimation  variance  is  given 

by  m2  =  20 .  When  >  7  ,  the  false  association  probability  is  seen  to  decrease  when  an  extra 

observable  state  is  used.  For  small  nu, ,  the  false  association  probability  increases  with  the 
number  of  the  observable  states.  However,  the  correct  association  probability  also  increases 
when  the  tracks  are  from  the  same  target  and  different  targets.  Figure  2(d)  shows  sgn(A^) 

versus  /(=  q2n)  and  /(=  m2n+x) .  Again  n  =  1.  When  the  number  of  observable  states  and  q2  are 
fixed,  it  can  be  seen  that  the  false  association  probability  decreases  when  the  ratio  m2n+x  of  the 
(/?+ 1  )th  state  difference  and  estimation  variance  is  above  certain  value. 
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3.  System  Model  for  Track  Association 


In  this  section,  we  present  a  new  modified  polar  coordinate  in  which  ESM/radar  track 
association  is  implemented.  The  observability  of  the  system  model  is  also  discussed. 

3.1  System  model 

The  geometry  of  the  new  MPC  is  shown  in  Figure  3.  The  new  MPC  is  a  combination  of 
Aidala’s  MPC  and  the  polar  coordinate  proposed  in  [18]  [19].  It  consists  of  the  5  states  that  are 
shown  to  be  observable  almost  anywhere  for  non-maneuvering  targets  and  observers.  The 
state  vector  of  the  Mi  target  is 

xk  Oo )  =  [0*  Oo )>  4 Oo ) /  rk  (4  )>  4  Oo )»  <4  Oo )» vk  (4 ) /  rk Oo )f >  k  =  1,  2,  ■■■,  K  . 

The  first  three  states  are  defined  in  Aidala’s  MPC  and  the  last  two  states  are  defined  in  the 
polar  coordinate  in  [18]  [19].  When  assuming  a  linearly  uniform  motion  model  for  the  target, 
the  trajectory  of  the  Mi  target  can  be  written  as  [18] 


**(0 

1 

g 

1 

-+>«» 

+ 

2 

i _ 

_%(0. 

_yM  +  0~t0)vky  _ 

(21) 


where  {xk{t),  yk{t)}  denotes  the  Cartesian  coordinates  of  the  target,  and  and  vky  are  the 
target  speeds  along  x-  and  y-axis,  respectively.  Aidala’s  MPC  can  be  written  as  [11] 

Qk  0)  =  °k  Oo )  +  tan-1  [s3  (t,t0)/s4(t,  t0 )] , 

•  /■,*  /  /vk  _  5i  0,  4  )ss  0>  0  +  s2  0, 4  K  0,  t0 ) 

rk  V l) '  rk  v l)  -  2  /  \  2  /  \ 

^3  ^0  )  +  54  ^0  ) 

•  _  ( t ,  t0  )s4  ( t ,  tQ)-s2  ( t ,  t0  )s3  ( t ,  t0 ) 

^k  \T )  ~  9  9  5 

53  O.^  +  ^OOo) 


(22) 

(23) 

(24) 


where 
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■'M (t, to)  =  0k(to)  +  [l/rk (t0 )] [ Wj (t,t0) cos dk (t0 ) -  w2 (t, t0 ) sin 6k (t0 )]  =  0(to )  , 
s2(t,to)  =  rk(to)/rk(to)  +  [l/rk(to)][wl(t,to)sm0k(to)  +  w2(t,to)cos0k(to)]  =  rk(to)/ rk (t0 ) , 
s3(t,to)  =  (t-to)&k(to)  +  [l/rk(to)][w3(t,to)cos0k(to)-w4(t,to)  sin  0k  (t0 )]  =  (t  - 10  )0(to ) , 

s4  (t,  t0  )  =  1  +  (t  - 10  )[rk  ( t0  )  /  rk  (70)]  (25) 


and 


s:0^wdA 

w,(t,t0) 

w2(t,t0) 

f  aky(A)dA 

w3(t,t0) 

\‘jt  -  A)kaJA)dA 

_w4(t,t0)_ 

\‘o(t  -  A)akv(A)clA 

(26) 


for  non-maneuvering  targets  and  observers,  a ^  and  aky  denote  the  relative  acceleration  of  the 
Mi  target  with  respect  to  the  ESM  sensor.  From  (22)  to  (26),  we  have 


^(0  =  6Wo)  +  tairl 


_ (i  t0  )0k  (tp ) _ 

l  +  (t-t0)[rk(t0)/rk(t0)] 


(27) 


Assume  that  t  =  t0+  ir  where  r  is  the  sensor  sampling  time  interval.  For  the  new  MPC,  the 
following  relationship  also  holds  [18] 


G  (0  =  tan”1  [  sin^-  +  /r  l>*  (*o  )/r*  Oo  )]  C0S[A-  Oo  )]  y 
k  1  cos[^  (t0 )]  +  ir  [vk  (t0  )/rk  (t0 )]  sin[^  (t0 )] 


(28) 


A(  0  -  A((o)  ’ 


vk( 0A*  0) 


_ ^OoV^Oo) _ 

4 1  +  2ir \yk  Oo  )  /  4  00 )]  sin[^  Oo  )  +  &  Oo  )]  +  0  r2  [v*  C^o  )  /  r<t  Oo  Xf 


(29) 

(30) 


We  model  the  ESM  sensor  bearing  measurement  as 
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Z*,«m(0  =  ^(0  +  «as»(0 


(31) 


where  nesm(t )  is  the  additive  measurement  noise  which  is  assumed  to  be  independently 
Gaussian  distributed  with  zero  mean.  The  radar  measurement  in  the  local  CC  is  given  by 


^k,  radar  (t)  =  h[xk(t)]  + 

ft radar  (0 


(32) 


where 

radar  (0  =  [**  (0,  4  (0.  yk  (0.  h  (Of , 
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ht[xk(t)] 

h2[xk{t)\_ 


0  [Xk  (0]  =  V**(0  +  .y*(0  =  yj(xk  (t0  )  +  irxk  (t0  ))2  +  ( yk  (t0  )  +  iryk  (/0  )  f 


h2  (xk  ( t ))  =  tan 


y*(fo)  +  frA(*o) 
**(*o)  + ”■**(*<>) 


(33) 


Based  on  (32),  we  can  estimate  the  target  tracks  in  the  CC  by  the  maximum  likelihood 
approach.  The  CC  to  the  new  MPC  conversion  can  be  done  by 


-l 
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rk(t)/rk(t)  = 

4(0  = 
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(34) 
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Note  that  the  states  of  the  proposed  MPC  are  not  redundant  because  no  state  can  be  computed 
from  the  rest. 

3.2  Observability  analysis 

We  use  the  following  rank  condition  matrix  to  test  the  observability  of  the  target  states 


o  = 


Ox 

0 


(35) 


where  O,  can  be  obtained  by  the  derivatives  of  0k  ( t )  in  Aidala’s  MPC  with  respect  to  its 
initial  state 


0\  (i,  y ) = i , 


Ot(2,j) 


r±  fa  xy  -  Qr  c°s(^-  (h + u  -  i)o  -  °k  (to )) 
rk  (to  +  (J  ~  (K) 


0\  (3,  j)  =  - 


rk  (t o  )0  ~  1)7  sin(6>A,  (t0  +  (y  -  l)r)  -  6>A,  (t0  )) 

h  +  O'  -  !)r) 


(36) 


and  02  can  be  obtained  from  the  derivatives  of  the  6k  (i t )  definition  in  (28), 


o2 


i 

0 

0 


d0k  (t0  +  r) 
d@k  (to) 

M (*o  +  T_) 
S(vk(t0)/rk(t0)) 

d&k  (to  +  0 


0^(to) 


dOk  (t0  +  2r) 

d0k(to) 
k_  (t0  +  2r) 
8(vk(to)/rk(to)) 
dOk  (t0  +  2t) 


dA(t0) 


(37) 


The  matrix  determinant  of  02  is  given  by  [18] 
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(38) 


-r  cos(<9a,  (f0  )  +  4  (ty  ))  -4r2  -  2r(vt  (f0  )  /  r±  (t0  ))  sin(flt  (f0  )  +  ^  (f0  )) 
‘ri  (*o )  ^2  (4  ) 

|  2rcos(flt(f0) +  &(*„))  -r2  - r(vt (*0 ) / rt (*0 )) sin(flt (f0 )  +  ^ (*0 )) 

4(0  ’  4(0 


where 


4  (h  )  =  C Vk  (to )  /  rk  (to  ))2  + 2-2  +  2r(vt  (*o )  /  rk  (to  ))  sin(^  (*o  )  +  <4  (*0  ))  > 

4  ('o )  =  O-  (4 )  1 1)-  (t0 ))  +  4r  +  4r(v/c  (t0 )  /  rk  (1() ))  sin(6^.  (70 )  +  (f)k  (t() )) .  (39) 

To  guarantee  det(02)  ^  0  ,  the  following  condition  must  be  satisfied, 

0*('o)  +  A('o)*f-  (40) 

To  guarantee  det(01 )  ^  0 ,  the  following  condition  must  be  satisfied, 

8k(t0+iT)-ek(t0)*Oorj.  (41) 

It  follows  that  the  states  in  the  new  MPC  are  observable  almost  anywhere  except  under  some 
very  special  conditions  expressed  as 


@k  (*o ) +  fik  (*o )  ~  2  ’ 


<4(4  +  O-4fc(4)  =  0ory  • 


(42) 


For  maneuvering  observers,  since  the  states  in  the  new  MPC  can  be  computed  from  the 
states  in  the  CC,  which  are  observable  [9],  they  are  observable. 
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4.  Maximum  Likelihood  Estimator 


In  this  section,  we  apply  the  ML  algorithm  for  estimating  the  tracks  in  the  new  MPC  for  the 
ESM  sensor.  The  CRBs  for  the  estimates  are  derived.  The  ESM/radar  association  logic  based 
on  the  state  estimates  and  their  respective  CRBs  is  also  discussed. 


4.1  ML  estimates  and  CRBs 


Because  the  new  MPC  is  a  combination  of  Aidala’s  MPC  and  the  polar  coordinate  in  [18], 
we  define  two  negative  log  likelihood  functions  separately.  The  first  is  based  on  the  ESM 
sensor  measurement  model  (1)  and  bearing  definition  (27)  in  Aidala’s  MPC 


Ll=N\og{nS2esm)  + 


N 

-Z 

i= 1 


o  +  it)  -  6k  (t0 )  -  tan 


ir&k(t0) 


l  +  ir[rk(t0)/rk(t0)\ 


.(43) 


where  the  constant  terms  are  ignored.  The  second  likelihood  is  based  on  the  ESM  sensor 
measurement  equation  (28)  and  bearing  definition  (25)  in  the  polar  coordinate  of  [18] 


l  n  2 

L2=N  log {n8lsm )  +  —  Z  (Zesm  (fo  + iT)  -  °k  (to  +  iT)) 

O esm 

=  N  log  {K8]sm )  +  -L  S  (zesm  (t0  +  iT)  (44) 

O esm  Z 

_tan-i  f  sint^  (f0 )]  +  iT  [ Vj_  (t0  )/rk  (f0 )]  cos[^  (f0 )]  |  ] 

1  cos[0t  (t0 )]  +  it  [vk  (t0 )/ rk  (t0 )]  sin [</>k  (t0 )]  J  y 

It  can  be  verified  that  minimizing  Lx  and  L2  is  equivalent  to  minimizing 


and 


A=  Z 


i= 1 


Zes,n  «0  +  iT)  ~  0k  (>0  )  ~  tai1  ' 


ir9k(t0) 


1\2 


l  +  ir[rk(t0)/rk(t0)] 


(45) 
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Ji  ='Z{zesm(t o  +ir)~  tan 

Z=1 


-1 


f  sin[^  fro  )]  +  K  fro  )/ r±  fro  )]  cos[4  fro  )]  1  ) 
cos[<9/(  (f0 )]  +  ir  [vk  fr0 )/ r*  (?„ )]  sin[^  (?0 )]  J  y 


(46) 


Obviously,  (45)  and  (46)  are  nonlinear  optimization  problems  and  numerical  techniques  are 
required.  The  Newton  optimization  method  is  used  to  estimate  the  Mi  target  state 

xk  fro )  =  [dk  fr0  ),  <t>k  fro  X 4  fro  X  vk  fro )  /  rk  fro  X  4  fro )  /  rk  fro  )f  •  The  (P+l)th  iteration  for 
calculating  Xk(t0)  is  given  by 


XiT  fro  )  =  X#  (to )  -  /u(p)HkjGkJ  (47) 

where  ,  (70  )  =  \Qk  fro  X  4  (t(j )  /  rk  (t() ),  6k  (t0  )]  ,  x k,  2  fro )  =  [  A-  fro  X  ft  fro  )  (  rk  fro  )]  >  /X  P '  is 
the  /?th  iteration  step-size,  <7^  t  is  the  gradient  vector  given  by 


N 


Gu  =2ZRkJ(i)rki(i) 


i= 1 


(48) 


where 


Yk,  1  (0  =  zesm  fro  + iT)  -  4  fro  )  -  tan  ' 


ir0k(to) 


rk,2  d)  =  f,™  (?0  +  iT)  ~  tan' 


!  +  /r[4(r0)/r)t(?o)] 

i  |  sin[gt  (A )]  +  ir  [vt  fr0 )/ rk  (f0 )]  cos[&  fr0 )] 

.  cos  [4  fro )]  +  ir  [vk  (to  )/4  (4 )]  sm[4  fro )] , 


(49) 


and 


**,,(0  = 


dfrt-fro  +*r) 

S-V./fro) 


(50) 


Using  the  bearing  definition  in  (27)  in  Aidala’s  MPC,  we  can  compute  7?^  (/)  by 
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d0kU()  +/r)  ^ 

S0*('o) 

86>A,(4+fr) _ (ir)240o) _ 

d(rk(t{i)/rk(t,)  [irOk  (t0 )]2  +  [1  + 1 irk (t0 ) / 4 (/0 )]2  ’ 

Q0  +  /r)  /r[l  +  izi\  (t0 ) /  r±  (t0 )]  ^ 

^a-(^o)  [it#*  (t0  )]2  +  [1  +  hrk  (f0  )  /  4  (t0  )]2 

/?,  2  (/)  can  be  calculated  based  on  the  bearing  definition  of  (28)  in  the  polar  coordinate  of 
[18]  as 


dflt  Oo  +  »r)  _ _ ~i2 ~  fo  Oo  )  /  4  Oo )]  sinlA  (4  )  +  <k  Oo )] _ 

Oo  )  [V*  (?o  )  /  4  O0  )]2  +  Or2  +  2i r[v/f  (t0  )  /  rk  (t0  )]  sin[<9A-  0(,  )  +  k  O0 )]  ’ 


d<400  +n~) _ -/rcos[4-(V)  + A-Oq)] _ 

9KOo)/^Oo)]  [v,00)/40o)]2  +4r2  +2ir[fA(t0)/4(t0)]sin[4(t0)  +  4(t0)] 


We  use  the  approximate  Hessian 


Hkl=j:Rkl(i)RTkl(i).  (53) 

1=1 

The  Fisher  information  matrix  (FIM)  of  the  ESM  track  estimates  for  the  Ath  target,  i.e., 
xk,\  Oo )  =  [0*  Oo  X  4  (4 ) 1  rk  Oo  X  4  (4  )f  and  Xt,2  Oo )  =  [  A  Oo  X  v/(  (t0 )  /  rk  (t0  )f  can  be  written 
as 


v  7  (=1 


of* 


k,l 


(/). 


(54) 


The  CRB  of  (70 )  is  the  diagonal  elements  of  the  inverse  FIM, 
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CRB, 


k,esm 


FIM 


k,  1 


FIMk\ 


(55) 


4.2  ML  estimator  and  CRB  for  radar 


The  negative  log  likelihood  function  of  the  radar  measurements  can  be  written  as 


where 


N 

£  =  -E  In 


+  \ll[zk,radar(to+iT)-h(xk(t0+iT)jQ-] 
^  i= 1 


27T^jdet(Q) 

[h, radar  (*0  +  *0  “  KXk (* 0  +  *0] 


(56) 


rv 

Ur,  radar  u 


J  6, radar 


The  target  state  x/(  radar(h)  can  be  estimated  similarly  by  using  the  Newton  optimization 

method.  When  the  radar  tracks  in  the  CC  are  estimated,  they  are  converted  into  the  new  MPC 
using  relationship  (34).  The  Fisher  Information  matrix  (FIM)  can  be  written  as 


FIM(Xk, radar 


dh{  (xk  (t0  +  ir)) 

T 

1 

0/z,  (vt(t0  +”-)) 

dxM 

O 

I _ 

“I 

tekVo) 

dh2(xk  (?0  +/r)) 

0  5l_ 

dh2(xk(t0  +  ir)) 

S^(t0) 

dxk(t0) 

(57) 


The  estimation  covariance  W  of  the  state  vector  xk  radar  (/,, )  =  [  xk  (/,, ),  xk  (t0 ),  yk  (/,, ),  yk  (/,,  )]r 
satisfies  the  following  inequality 


W  -E[xkmdar  Xk,radar][*k, radar  Xk,radaA  -  FIM  (xkmdar).  (58) 
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The  FIM  of  Xk  (t)  in  the  new  MPC  can  be  derived  from  the  FIM  of  radar  track  in  the  CC 


as  [10] 


FIM radar  (Xk  (t0  )) 


N 

I 


i=l 


I T 


dX k_  (4  +  it) 

^Xk,i radar  ({0  ) 


dfh(xk(t0+iT)) 

dxk(t0) 

dhi  (Xk  (*o  +  /r  » 
dxk(t0) 


1 

3/i,  (xk  (t0  +  it)) 

~s2r  0 " 

—  1 

dxk(t0) 

dXk  (t0  +  it) 

.0  sl_ 

dh2(xk(t0  +/'r)) 
dxk(t0) 

^Xk, radar  (?0  ) 

(59) 


The  CRB  for  Xk  (t0 )  is  given  by 


CRBk,adar  =[FIMmdar(Xk(t0)] 


(60) 


4.3  ESM/radar  track  association 


Consider  two  tracks  Xradar  and  Xesm  by  the  ML  estimators.  Under  the  hypothesis  that 
Xradar  and  Xesm  are  independent  and  Gaussian  distributed,  the  association  decision  test  is 
defined  by 


( ^ radar  ~  ^esm  )  (  CRB  radar  +  CRBesm  )  ( ^ radar  ~  ^esm  )  ~T  (6 1 ) 

where  T  is  a  preset  threshold.  The  radar  and  ESM  tracks  are  said  to  be  associated  when  (61) 
holds.  Otherwise,  they  are  considered  to  be  from  different  targets.  In  practice,  since  the  real 
target  states  are  not  known,  we  use  the  estimated  target  states  and  equations  (54)  and  (60)  to 
evaluate  the  CRB  for  the  target  state  estimates  in  (61). 
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5.  Experimental  Results 


In  this  section,  we  use  computer  simulations  to  assess  the  ESM/radar  track  association 
performance  in  the  new  MPC  for  different  scenarios.  These  scenarios  include:  (1)  a  single 
target  with  locally  linear  dynamic  model,  (2)  two  close  targets  with  locally  linear  dynamic 
models,  which  are  moving  away  from  the  radar  and  ESM  sensor,  (3)  two  close  targets  with 
locally  linear  dynamic  models,  which  are  moving  towards  the  sensors,  (4)  two  close  targets 
with  locally  linear  dynamic  models,  which  are  moving  away  from  the  sensors,  and  (5)  two 
targets  with  locally  linear  dynamic  model,  and  the  observers  are  maneuvering.  In  scenarios  (1) 
to  (4),  the  observers  are  assumed  to  be  non-maneuvering. 

Figure  4  compares  the  variances  of  the  ML  estimates  of  the  ESM  track  with  the  CRBs.  The 
initial  position  of  the  target  is  given  by  (-8000  m,  12000  m).  The  initial  velocities  along  x-  and 
y-axis  are  320  mis  and  230  m/s ,  respectively.  Figures  4(a)  to  4(e)  show  the  estimation  errors 
of  the  five  components  of  the  state  vector,  respectively.  It  can  be  seen  that  the  variances  of  the 
ML  estimates  of  the  ESM  tracks  approach  their  CRBs  as  the  number  of  measurements 
increases. 

Figure  5  and  Figure  6  show  the  performance  of  different  association  algorithms  using  the 
proposed  MPC  with  different  testing  thresholds.  In  Figure  5,  the  testing  threshold  T  is  set  to 
3.84  (corresponding  to  the  correct  association  probability  pc  =  0.95  for  a  chi-square 
distribution  of  freedom  1).  Three  association  algorithms  are  compared.  The  first  uses  only  one 
state  (bearing)  as  in  [6]  [8].  The  second  approach  uses  two  states  (bearing  and  range 
rate/range)  and  the  third  one  combines  all  the  five  observable  states  in  the  new  MPC.  In 
Figure  5(a),  the  testing  variable  for  a  given  number  of  samples,  k ,  is  defined  as 

y  =  [Ktar  Km  ){CRBradar  (k)  +  CRBsm  {k))^  (Xradar  Km  )'  • 


In  Figure  5(a),  the  correct  association  probabilities  with  one,  two,  and  five  states  are  about 
0.95,  0.98,  and  0.995,  respectively.  The  performance  of  the  ESM/radar  track  association 
indeed  improves  with  the  use  of  more  observable  states.  It  is  also  observed  that,  when  the 
number  of  measurements  increases,  both  the  target  state  estimation  errors  and  the  CRB 
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decrease  while  the  correct  association  probability  keeps  almost  the  same.  In  Figure  5(b),  the 
testing  variable  for  k  is  defined  as 


y=(x„^-xJj{cRB^m+cRB^m)-1  (xmdar-Km)T  ■ 

It  can  be  seen  that  the  correct  association  probability  improves  with  the  increase  of  the 
number  of  measurements  because  of  the  smaller  estimation  errors  of  {Xradar  -X^j.  In 

Figure  6,  T  is  set  to  2.05  (corresponding  to  a  correct  association  probability  pc  =0.85  for  a 
chi-square  distribution  of  freedom  1). 

Figure  7  shows  the  correct  association  probabilities  of  different  association  techniques 
when  two  targets  are  moving  away  from  the  observer.  Both  targets  are  observed  by  the  radar 
and  ESM  sensor.  The  initial  states  of  target  A  and  target  B  are  (-8000  m,  320  m/s ,  12000  m, 
230  m/s)  and  (-8000  m,  230  m/s ,  10000  m,  320  m/s),  respectively.  T  set  to  be  3.84.  Let  pcaa 
and  pcbb  denote  the  probabilities  of  correctly  associating  ESM  track  A  with  radar  track  A, 
and  ESM  track  B  with  radar  track  B,  respectively,  and  pcab  and  pcba  represent  the 
probabilities  of  not  associating  ESM  track  A  with  radar  track  B,  ESM  track  B  with  radar  track 
A,  respectively.  The  correct  association  probability  in  Figure  7(a)  and  7(b)  is  defined  as 
Pc  =  ( Pcaa  +  Pcbb  +  Pcab  +  Pcba )  /  4  .  The  definitions  of  the  testing  variable  Y  in  Figure  7(a)  and 
7(b)  are  the  same  as  in  Figure  5(a)  and  5(b),  respectively.  Figures  7  (a)  and  7(b)  show  that  the 
probability  of  correctly  associating  both  the  same  target  and  different  targets  can  be  improved 
effectively  by  using  more  observable  states.  Another  index  for  evaluating  the  performance  of 
association  is  the  probability  of  correctly  associating  the  same  target,  which  is  defined  as 
Pc  =(Pcaa  +  Pcbb) /2,  and  shown  in  Figure  7(c)  and  7(d).  From  Figure  7(c)  and  7(d),  we  can 
observe  than  the  probability  of  correctly  associating  the  same  target  increases  when  more 
observable  states  are  used. 

Figure  8  shows  the  probabilities  of  falsely  associating  one  target  from  another.  The 
simulation  parameters  are  the  same  as  in  Figure  7.  The  false  association  probability  in  Figure 
8  is  defined  as  pf  =  (p  fah  +  pfha )  /  2  ,  where  pfab  and  p denote  the  probabilities  of  falsely 

associating  ESM  track  A  with  radar  track  B,  and  ESM  track  B  with  radar  Track  A, 
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respectively.  Figure  8  shows  that  the  false  association  probabilities  based  on  two  states 
(bearing  and  heading  direction)  and  five  states  are  lower  than  that  of  the  method  based  on 
bearing  only.  In  Figures  7  and  8,  it  is  shown  that  the  two  states  and  the  five  states  based 
methods  not  only  increase  the  probabilities  of  correctly  associating  tracks  from  the  same 
target,  but  also  decrease  the  probabilities  of  falsely  associating  tracks  from  different  targets  at 
the  same  time. 

Figure  9  shows  the  correct  association  probabilities  of  different  association  techniques 
when  two  targets  are  moving  towards  the  observers.  The  initial  states  of  target  A  and  B  are 
(20000  m ,  100  m/s ,  30000  m,  200  m/s)  and  (20000  m,  150  m/s ,  25000  m,  200  m/s), 
respectively.  T  is  set  to  be  3.84.  The  correct  association  probability  in  Figure  9(a)  and  9(b)  is 
defined  as  pc  =  (pcaa  +  pcbb  +  pcab  +  pcba)/ 4  .  The  test  Y  in  Figure  9(a)  and  9(b)  is  the  same 
as  in  Figure  5(a)  and  5(b),  respectively.  In  Figures  9  (a)  and  9(b),  it  is  shown  that  the 
probability  of  correctly  associating  both  the  same  target  and  different  targets  improves 
effectively  by  using  more  observable  states.  The  correct  association  probability  in  Figures 
9(c)  and  9(d)  is  defined  as  pc  =  (pcaa  +  pcbh)/2.  Figure  10  shows  the  probabilities  of  falsely 
associating  one  target  from  another.  The  simulation  parameters  are  the  same  as  in  Figure  9.  In 
this  scenario,  the  false  association  probability  decreases  when  more  observable  states  are 
used. 

Figure  11  shows  the  correct  association  probabilities  for  different  association  techniques 
when  two  close  targets  are  moving  away  from  the  observers.  The  initial  states  of  Target  A  and 
Target  B  are  (-8000  m ,  320  m/s ,  12000  m ,  230  m/s)  and  (-8000  m,  320  m/s,  11000  m,  230 
m/s),  respectively.  T  is  3.84  for  Figure  11(a)  and  11(b),  and  is  2.05  for  Figures  11(c)  and 
11(d).  In  Figure  11(a)  and  7(c),  the  correct  association  probability  is  given  by 
Pc  =(Pcaa  +  Pcbb  +  Pcab  +  Pcba)1 4  •  The  test  Y  in  Figure  11(a)  and  11(b)  is  the  same  as  in 
Figure  5(a)  and  5(b),  respectively.  Figures  11(b)  and  11(d)  show  that  the  probability  of 
correctly  associating  the  same  target  improves  in  the  case  of  two  close  targets  moving  away 
from  the  observers.  In  Figures  1 1(a)  and  1 1(b),  it  is  observed  that  the  total  correct  association 
probability  pc  =  (pcaa  +  pcbb  +  pcab  +  pcba )  /  4  decreases  for  small  numbers  of  measurements 
when  more  observable  states  are  used.  However,  it  improves  as  the  number  of  measurements 
increases.  The  five  states  based  association  becomes  the  best  when  the  number  of 
measurements  is  70  for  T=3.84,  and  40  for  the  T=2.05. 


23 


Figure  12  shows  the  probabilities  of  falsely  associating  one  target  from  another.  The 
simulation  parameters  are  the  same  as  in  Figure  1 1.  In  Figure  12(a)  and  12(b),  T is  set  to  3.84 
and  2.05,  respectively.  When  the  number  of  measurements  is  small,  it  is  observed  that  the 
false  association  probability  increases  with  the  use  of  extra  observable  states,  which  causes 
the  deterioration  of  the  total  correct  association  probability  shown  in  Figures  1 1(a)  and  1 1(c). 
The  reason  is  that  the  ratio  of  the  state  difference  and  estimation  variance  is  small  with  small 
number  of  measurements.  When  the  number  of  measurements  increases,  the  probability  of 
falsely  associating  two  close  targets  approaches  0,  and  the  total  association  probability  with 
more  states  begin  to  its  advantages.  Although  the  false  association  probability  may  increases 
on  its  own,  the  probability  of  correctly  associating  the  same  target  always  improves  with  the 
use  of  extra  observable  states,  as  shown  in  Figures  1 1(b)  and  1 1(d). 

Figure  13  shows  the  correct  association  probabilities  for  maneuvering  observers.  The 
tracks  of  the  targets  and  observers  are  shown  in  Figure  13(a).  The  correct  association 
probabilities  for  different  association  techniques  are  shown  in  Figure  13(b).  When  the 
observer  is  maneuvering,  the  states  of  the  targets  in  CC,  i.e.,  positions  and  velocities  along  x- 
and  y-axis  are  observable.  We  use  an  ML  estimator,  which  is  similar  to  that  introduced  in 
Section  5,  for  estimating  the  radar  tracks  in  the  CC.  The  estimated  tracks  are  then  converted  to 
the  new  MPC  for  association.  Figure  13(b)  shows  that  the  correct  association  probability 
improves  with  the  use  of  extra  observable  states  for  maneuvering  observers. 
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6.  Conclusions 


In  this  report,  a  new  modified  polar  coordinate  (MPC)  system  has  been  proposed  based  on  a 
theoretical  performance  analysis  of  the  classical  track-to-track  association  techniques  for 
maneuvering  and  non-maneuvering  targets  and  observers.  It  is  shown  that  the  correct 
association  probabilities  will  decrease  with  the  introduction  of  bias  into  the  observable  states 
non-maneuvering  targets  and  observers.  The  performance  of  the  ESM/radar  track  association 
can  be  improved  effectively  if  more  observable  states  are  used  for  association.  Based  on  this 
conclusion,  we  proposed  ESM/radar  track  association  algorithm  in  the  new  MPC  which 
includes  extra  observable  states  in  addition  to  Aidala’s  MPC.  A  maximum  likelihood  (ML) 
algorithm  is  proposed  to  estimate  the  ESM  tracks  in  the  new  MPC.  The  new  MPC  formulation 
is  applicable  to  both  maneuvering  and  non-maneuvering  targets  and  observers.  The  CRBs  of 
the  radar  and  ESM  track  estimates  are  also  derived.  Computer  simulations  were  used  to 
demonstrate  the  performance  of  the  proposed  algorithm.  It  was  shown  that  the  track 
association  performance  improved  effectively  by  using  more  observable  states  in  the  new 
MPC  in  different  scenarios. 
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Figure  1.  Figure  1.  Correct  association  probabilities  under  different  cases:  (a)  The 
radar  and  ESM  tracks  are  related  to  the  same  target;  (b)  The  radar  and  ESM  tracks 
are  related  to  different  targets.  The  pc  atw  =  0  in  (a)  and  (b)  shows  the  correct 
association  probability  in  case  1  and  case  3,  respectively,  i.e.  no  estimation  biases  for 

the  target  states. 
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(a)  (b) 


Figure  2.  Association  probability  versus  the  number  of  observable  states  used  for 
ESM/radar  track  association.  The  testing  threshold  is  fixed  to  be  T  =  3.84.  (a) 
probability  of  correctly  associating  tracks  from  the  same  target;  (b)  probability  of 
falsely  associating  tracks  from  different  targets;  (c)  probability  of  correctly  associating 

tracks  from  the  same  targets  and  from  different  targets;  (d)  sgn(A pf )  versus  ?'(=  q2n ) 

and  j(=m2n+x) ,  n  =  1  here. 
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(a) _  _ (b) 


(c)  (d) 


(e) 


Figure  4.  Comparison  of  the  ML  estimator  and  CRB  for  ESM  tracking:  (a)  bearing;  (b) 
range  rate/range;  (c)  bearing  rate;  (d)  heading  direction;  (e)  target  velocity/range. 
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Figure  5.  Correct  association  probability  for  the  single  target.  The  testing  threshold  is 
3.84.  (a)  The  variance  in  the  testing  variable  changes  with  the  sample  number;  (b) 
The  variance  in  the  testing  variable  is  fixed. 


Figure  6.  Correct  association  probability  for  the  single  target.  The  testing  threshold  is 
2.05.  (a)  The  variance  in  the  testing  variable  changes  with  the  sample  number;  (b) 
The  variance  in  the  testing  variable  is  fixed. 
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Figure  7.  Correct  association  probabilities  for  two  close  targets  moving  away  from  the 
observer:  (a)  Pc  =  ( Pcaa  +  Pcbb  +  Pcab  +  Pcba)/4.  Changed  variance  in  the  testing 
variable;  (b)  Pc  =  ( Pcaa  +  Pchh )  12;  Fixed  variance  in  the  testing  variable;  (c) 

Pc  —  ^ Pcaa  +  Pcbb  +  Pcab  +  Pcba) 1 4  ■  Changed  variance  in  the  testing  variable;  (d) 

Pc  =  ( +  PCbb )  / 2  ■  Fixed  variance  in  the  testing  variable. 
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Figure  8.  False  association  probabilities  for  the  two  targets  in  Figure  7. 
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(a) 


(b) 


(b)  (d) 

Figure  9.  Correct  association  probabilities  for  two  close  targets  moving  toward  the 
observer:  (a)  pc  =  ( pcaa  +  pcbb  +  pcab  +  pcba)i 4.  Changed  variance  in  the  testing 

variable;  (b)  pc  =  {pcaa  +  pcbb)/2;  Fixed  variance  in  the  testing  variable;  (c) 

Pc  =  ( Pcaa  +  Pcbb  +  Pcab  +  Pcba) ' 4  ■  Changed  variance  in  the  testing  variable;  (d) 

Pc  =  ( Pcaa  +  Pcbb )  / 2  ■  Fixed  variance  in  the  testing  variable. 
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Figure  10.  False  association  probabilities  for  the  two  targets  in  Figure  9. 
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(a)  (b) 


(c) 


(d) 


Figure  1 1.  Correct  association  probabilities  for  two  very  close  targets  moving  away 
from  the  observer:  (a)  T  =  3.84,  Pc  =  (Pcaa  +  Pcbb  +  Pcab  +  Pcba)/ 4;  (b)  T  =  2.05, 

Pc  =(Pcaa  +  Pcbb)12- 
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Figure  12.  False  association  probabilities  for  the  two  targets  in  Figure  1 1:  (a)  T  = 

3.84;  (b)  T  =  2.05. 


(a)  (b) 

Figure  13.  Correct  association  probabilities  for  maneuvering  observer:  (a)  tracks  for 
targets  and  observers;  (b)  correct  association  probabilities. 

Pc  =  ( Pcoa  +  Pcbb  +  Pcab  +  Pcba) ^  ■  The  testin9  threshold  is  2.05. 
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